Response of upper tropospheric water vapor to global warming and ENSO

The upper tropospheric water vapor is a key component of Earth's climate. Understanding variations in upper tropospheric water vapor and identifying its influencing factors is crucial for enhancing our comprehension of global climate change. While many studies have shown the impact of El Niño-Southern Oscillation (ENSO) and global warming on water vapor, how they affect the upper tropospheric water vapor remains unclear. Long-term, high-precision ERA5 specific humidity data from the European Centre for Medium-Range Weather Forecasts (ECMWF) provided the data foundation for this study. On this basis, we successfully obtained the patterns of global warming (Independent Component 1, IC1) and ENSO (Independent Component 2, IC2) by employing the strategy of independent component analysis (ICA) combined with non-parametric optimal dimension selection to investigate the upper tropospheric water vapor variations and responses to ENSO and global warming. The results indicate that global warming and ENSO are the primary factors contributing to water vapor variations in the upper troposphere, achieving the significant correlations of 0.87 and 0.61 with water vapor anomalies respectively. Together, they account for 86% of the global interannual variations in water vapor. Consistent with previous studies, our findings also find positive anomalies in upper tropospheric water vapor during El Niño years and negative anomalies during La Niña years. Moreover, the influence extent of ENSO on upper tropospheric water vapor varies with the changing seasons.


Data and methods
The accuracy of ERA5 water vapor data relies on data assimilation methods and the stability and precision of observational systems.Microwave observations and infrared datasets (e.g., AMSU, HIRS) provide accurate estimates of water vapor variations in the upper troposphere 12,31,32 , while spectral or radio occultation data (e.g., AIRS) provides information on the distribution of specific humidity and relative humidity 33,34 .Since satellite data are important for determining water vapor changes, especially over the oceans, the scarcity of satellite data prior to the mid-1990s resulted in unrealistic trends in the ERA5 dataset 12,35 .The specific humidity data used in this article stems from the fifth generation of the ECMWF atmospheric reanalysis data, which is available in the European Union's Earth observation program ERA5 Copernicus Climate Data Store (https:// cds.clima te.coper nicus.eu).ERA 5 specific humidity dataset is a product of combining data assimilation with weather forecast modelling system with high resolution and accuracy, which is crucial for studying the global water cycle, climate change, and extreme weather events 36 .The high spatial resolution (0.25° × 0.25°) and high temporal resolution (hourly) of ERA5 data make it possible to conduct precise studies of trends and variations in atmospheric water vapor 15 .A long-term time series of ERA5 specific humidity data spanning from January 1959 to December 2021 is selected in this paper to investigate the effects of ENSO and global warming on atmospheric water vapor in both temporal and spatial domains.
ENSO, driven by sea surface temperature variations in the central and eastern equatorial Pacific, is an important source of tropical Pacific climate anomalies.It not only has a significant impact on water vapor, but also is one of the most important factors contributing to global climate change.The Oceanic Niño Index (ONI) is the most suitable indicator for characterizing ENSO events 37 , and it is also one of the most commonly used indices for defining El Niño and La Niña events.Therefore, in order to quantitatively analyze the temporal and spatial effects of ENSO on water vapor, ONI (https:// psl.noaa.gov/ data/ corre lation/ oni.data) from the National Oceanic and Atmospheric Administration (NOAA) Physical Sciences Laboratory is selected in this paper.
Global warming is a long-term phenomenon characterized by the widespread increase in the Earth's surface temperature caused by greenhouse gases.It has profound impacts on ecosystems, economies, and societies 38,39 .The global mean surface temperature (GMST) is one of the most important indicators for measuring global warming and climate change 40 .Therefore, GMST (https:// clima te.nasa.gov/ vital-signs/ global-tempe rature/) is selected from (National Aeronautics and Space Administration) NASA Goddard Institute for Space Studies (GISS) to investigate the temporal and spatial response of water vapor to global warming.
Climate time series exhibit nonlinear characteristics.In order to overcome the impact of nonlinearity on the quality of component extraction, a new processing strategy of the optimal lowpass filtering combined with ICA is proposed for the first time.The complete process of the proposed method is illustrated in Fig. 1.
The process flow can be divided into two parts: pretreatment and component extraction.The pretreatment primarily consists of the following four steps.Firstly, the grid specific humidity data is obtained from the ERA5 dataset.Secondly, the specific humidity anomalies signals were obtained by eliminating the annual and monthto-month variations 41 .In detail, the annual cycle was eliminated by subtracting the month-by-month climatic mean state.On the other hand, the 1-2-1 three point binomial filter was used to eliminate the monthly variability 42 .Thirdly, the optimal output dimension is determined by the non-parametric method based on the specific humidity anomalies signals 30,43 .The optimal output dimension is 2 according to the data used in this paper.Fourthly, in order to eliminate the noise from specific humidity anomalies signals, the optimal lowpass filtering was applied 44 .However, due to the uncertainty of the optimal cut-off frequency of lowpass filtering, it is essential to compare the results with different cut-off frequencies to determine the optimal cut-off frequency.
According to the selected optimal output dimension, independent components can be estimated by ICA for all filtering cut-off frequencies.For each filtering cut-off frequency, Pearson correlation analysis was carried out between the ICs and ONI/GMST.The optimal cut-off frequency of the lowpass filter was determined when the maximum absolute value of the correlation coefficient was obtained.The time series of ICs which have the www.nature.com/scientificreports/maximum absolute value of the correlation coefficient with ONI were regarded as ENSO pattern.Similarly, the time series of ICs which have the maximum absolute value of the correlation coefficient with GMST were regarded as global warming pattern.

Characteristics of water vapor variations and the affecting factors
To provide a clearer analysis of the characteristics of upper tropospheric water vapor variations and assess the factors influencing water vapor variations, the results of the time series after each processing step of the component extraction flow chart are presented here.Figure 2a displays the global mean specific humidity data, indicating a noticeable seasonal cycle in the specific humidity data.This is consistent with the findings of Patel et al. 17 , whose research investigated the long-term variations in global tropospheric water vapor and identified significant seasonal variations.This study also found that water vapor content is highest in summer and lowest in winter 17 .In addition, to investigate the factors influencing water vapor variations, this study used the strategy of ICA combined with the non-parametric optimal dimension selection to separate the corresponding patterns of global warming and ENSO from specific humidity anomalies.The results are shown in Fig. 2e and f, respectively.
From the time series of specific humidity anomalies after processing (Fig. 2b, c), it is evident that the time series exhibits a slow upward trend, which is consistent with Dessler et al. 1 , who also found a long-term upward trend in the upper tropospheric water vapor data.This may be associated with the increase in deep convection caused by global temperature rise 17 .From Fig. 2b and c, it can be observed that the majority of strong El Niño years correspond to the maximum values in the specific humidity anomalies time series, while strong La Niña years correspond to minimum values.This is consistent with the findings of Saha et al. 4 and Tian et al. 45 , where the peak values of water vapor anomalies are associated with El Niño years, while the minimum values are www.nature.com/scientificreports/associated with La Niña years.There may be a lag in the response of upper tropospheric water vapor to ENSO phenomenon, causing some ENSO years not having a good one-to-one correspondence with the maximum or minimum values.The lag phenomenon can be attributed to the fact that sea surface temperatures reach their peak in winter (DJF = December, January, and February) during ENSO events.Sea surface temperatures alters the pressure gradients over the Pacific Ocean, affecting trade winds and the Walker circulation.Subsequently, the Walker circulation determines the position of deep convection in the Pacific, causing the Kelvin wave to take 2-3 months to cross the Pacific.Additionally, the transition of convection also takes time.Only after the convection changes can water vapor be transported to the upper troposphere.Hence, there is a lag in the impact of ENSO on upper tropospheric water vapor, requiring several months to respond to ENSO 4 .The spatial distribution of water vapor flux exhibits clear latitude characteristics, with water vapor flux between oceans and land decreasing as latitude increases 18 .Additionally, under the influence of global warming, the spatial extent and the feedback of water vapor becomes stronger over time 23 .Figure 3a shows the latitudetime evolution of global specific humidity.It is evident that there is a significant annual cycle, and the influence range and intensity of water vapor increases over time.This result is consistent with the findings in Fig. 2a, which are associated with the positive feedback mechanisms between global warming and water vapor 17 .Furthermore, it can also be seen in Fig. 3a that upper tropospheric water vapor is primarily concentrated in tropical regions (30° S-30° N), which is associated with the fact that the tropics contribute to approximately 2/3 of the water vapor feedback in the upper troposphere 6 .Figures 3b represent the latitude-time evolution of specific humidity anomalies after removing the annual cycle and monthly variations, respectively.It can be observed that the water vapor anomalies were primarily negative before the twenty-first century, while a positive response became more dominant after the twenty-first century.This may be attributed to the upward trend in water vapor anomalies under the influence of global warming, which results in specific humidity data before the twenty-first century generally remained below average, while data after the twenty-first century are almost always above average.Figure 3c  The study of Shine et al. 46 points out that from a radiometric perspective, percentage variations in water vapor rather than absolute variations are more relevant to water vapor feedback.The findings of our study indicate that 53% of the global interannual variations in water vapor can be attributed to global warming, while 33% of the global interannual variations in water vapor can be explained by ENSO.In addition, the remaining 14% of the water vapor change may be related to water vapor transport caused by wind, seasonal changes, and other factors 47 .Moreover, to further investigate the latitude characteristics of water vapor and analyze the contributions of global warming and ENSO to water vapor, the latitude distribution of the deseasonalized water vapor anomalies, global warming pattern (IC1), ENSO pattern (IC2), and global warming plus ENSO pattern were generated (Fig. 4).
From the results, it can be observed that the water vapor anomalies associated with global warming pattern (IC1) largely capture variations of water vapor in the upper troposphere.This is consistent with the fact that the IC1 pattern explains 53% of the global interannual variations in water vapor, indicating that global warming pattern is a primary factor driving variations in water vapor in the upper troposphere.Furthermore, we observed that the influence of ENSO pattern (IC2) on water vapor anomalies in the upper troposphere is primarily concentrated in tropical regions.This may be linked to the fact that ENSO originates in the tropical Pacific and tropical convection 4,48 .Moreover, we combined the contributions of global warming pattern (IC1) and ENSO pattern (IC2) on water vapor in the upper troposphere.We found that, except for a few areas (e.g., 30° S-60° S), the combined results are closely consistent with the upper tropospheric water vapor anomalies.The changes in water vapor result from the influence of multiple factors, while IC1 and IC2 respectively representing global warming and ENSO only capture the impacts of these two factors on water vapor.Besides, global warming is a steady, globally uniform forcing phenomenon, whereas the impact of ENSO on water vapor is primarily existed in tropical regions.Thus, the mentioned reasons could result in differences in the relationship between the deseasonalized water vapor anomalies and IC1, IC2, and IC1 + IC2 in certain regions.Therefore, we believe that global warming and ENSO patterns are significant factors contributing to variations in the upper tropospheric water vapor.

The response of the upper tropospheric water vapor to global warming
To quantitatively analyze the relationship between upper tropospheric water vapor variations and global warming, the time series of the first independent component (IC1) of upper tropospheric specific humidity anomalies is obtained by employing the strategy of ICA combined with the non-parametric optimal dimension selection, as shown in Fig. 5.
The results of Independent Component Analysis indicate that the pattern of IC1 accounts for 53% of the global interannual variations in water vapor.The correlation coefficient between the time series of IC1 and the upper tropospheric water vapor anomalies is 0.87 (Fig. 5a).At the same time, it exhibits a significant correlation with GMST indicator representing the global warming signal, with a correlation coefficient as high as 0.91 (Fig. 5b) with a lag of 11 months (i.e., the maximum correlation coefficient with the GMST indicator is achieved after 11 months).As global warming is the response to a steady, globally uniform forcing 23 , altering convection is slower.Therefore, it requires a longer time to impact upper tropospheric water vapor compared to ENSO.The above results indicate that the pattern of IC1 captures the main features of water vapor variations, and a primary factor leading to upper tropospheric water vapor variations is global warming.Furthermore, we observed that there is usually a local peak or minimum in both water vapor anomalies and GMST during strong El Niño or La Niña years (highlighted in red and blue backgrounds in the figures).This suggests that strong ENSO events may have a significant impact on both global warming and water vapor anomalies.A possible explanation is that ENSO alters tropical convection, leading to variations in upper tropospheric water vapor.Subsequently, under the influence of the positive feedback mechanism of water vapor, it contributes to the variations in greenhouse gases, thereby influencing global warming.The correlation coefficient between the global warming pattern and the upper tropospheric water vapor anomalies is 0.87.Moreover, the global warming pattern also explains as much as 53% of the global interannual variations in water vapor, indicating that the upper tropospheric water vapor change is highly dependent on global warming.This correlation result also supports the results of Allan et al. (2022) that at higher altitudes in the troposphere, the global mean water vapor responds more strongly to temperature changes.This may be related to the following reason.According to the Clausius-Clapeyron equation, atmospheric water vapor increases with rising temperatures 49 .Therefore, the atmospheric water vapor content will increase under the influence of global warming.Water vapor is also a greenhouse gas, and with the increase in greenhouse gases, the amount of longwave radiation reaching the Earth's surface also intensifies, further exacerbating global warming 23 .Consequently, this leads to a heightened correlation between global warming patterns and water vapor in the upper troposphere.www.nature.com/scientificreports/Global warming intensifies the water vapor transport between the ocean and land 18 .To analyze the spatial response of upper tropospheric water vapor to global warming, Fig. 6 depicts the global distribution of the 63-year long-term mean of upper tropospheric specific humidity (Fig. 6a) and the spatial pattern of IC1 (Fig. 6b).The results of the global distribution of long-term mean of upper tropospheric specific humidity (Fig. 6a) are consistent with the water vapor distribution obtained by Takahashi et al. 50and Tian et al. 51 .The distribution of upper tropospheric water vapor shown in Fig. 6a is very close to the spatial pattern of IC1 (global warming pattern), indicating that the upper tropospheric water vapor is significantly influenced by the global warming pattern.Wang et al. 18 pointed out that the exchange of water vapor between oceans and land near South America, Africa, and Australia is more sensitive under global warming conditions, and the spatial pattern of IC1 (Fig. 6b) also illustrates this spatial characteristic.The Intergovernmental Panel on Climate Change (IPCC) has indicated  that global warming will further strengthen the large-scale atmospheric circulation, leading to significant changes in the global and regional water cycles 52 .This will intensify the transport of water vapor between oceans and land, consequently increasing the likelihood of extreme precipitation events and floods 18 .To investigate the relationship between global warming patterns and water vapor transport, this study employs the water vapor flux to depict the exchange of water vapor between the oceans and the land (Fig. 6c).The results indicate that the region of higher water vapor flux is approximately the same as the region affected by global warming (60° N-60° S).Additionally, regions such as Southeast Asia, and the lower part of North America experience pronounced warming (Fig. 6b), while at the same time, water vapor fluxes are also significant where land borders the ocean near these locations.Simultaneously, there is a significant water vapor flux in the mid-latitude region of the southern hemisphere pacific (30° S-60° S), which may be associated with the expansion of the Hadley circulation in the southern hemisphere, influenced by global warming (IPCC).Above results indirectly validate that under the influence of global warming significant changes in the global and regional hydrological cycles are occurring, and that water vapor exchange between oceans and land will also be intensified 18,52,53 .

The response of the upper tropospheric water vapor to ENSO
ENSO exhibits a strong influence on the interannual anomalies of tropical upper tropospheric water vapor 4 .Therefore, to quantitatively analyze the response of upper tropospheric water vapor to ENSO, the time series of another independent component (IC2) of upper tropospheric specific humidity anomalies is also obtained, as shown in Fig. 7.
According to the results of independent component analysis, the pattern of IC2 accounts for 33% of the global interannual variations in water vapor.Due to the influence of the linear trend caused by global warming on upper tropospheric water vapor anomalies, the upper tropospheric water vapor anomalies were detrended according to the method of Lu et al. 16 .The correlation coefficient between the time series of upper tropospheric water vapor anomalies after removing the linear trend and IC2 is 0.61 (Fig. 7a).Additionally, with an 8 month lag, a significant correlation was found between ONI and a correlation coefficient as high as 0.91 (Fig. 7b).The above results indicate that the pattern of IC2 also captures the characteristics of water vapor variations and that the other major factor contributing to water vapor variations in the upper troposphere is ENSO.
ENSO usually peaks during winter, which is also the period when ENSO has the strongest impact on water vapor 4 .Therefore, to analyze the spatial response of upper tropospheric water vapor to ENSO, Fig. 8 illustrates the El Niño-year mean of upper tropospheric specific humidity during the DJF months (Fig. 8a), the La Niñayear mean of upper tropospheric specific humidity at the DJF months (Fig. 8b), and the spatial pattern of IC2 (Fig. 8c).The spatial pattern results indicate that the global distribution of upper tropospheric specific humidity is mainly positive response during El Niño events, while the global distribution of upper tropospheric specific humidity is significant negative response during La Niña events.This result is consistent with the results of Saha et al. 4 .Additionally, we found that the El Niño-year mean of upper tropospheric specific humidity during the DJF months (Fig. 8a) is very consistent with the spatial pattern of IC2, and both of them exhibit the classic ENSO pattern 48 , indicating that the upper tropospheric water vapor is affected by the El Niño pattern during El Niño events, resulting in a consistent spatial response.Generally, ENSO has a positive response on upper tropospheric water vapor under the influence of El Niño, leading to relatively moist conditions.Similarly, ENSO has a negative response on upper tropospheric water vapor during La Niña years, resulting in relatively dry conditions in the upper troposphere.In addition, relatively dry patches (blue patches) appeared next to warm pool under El Niño conditions (Fig. 8a, c) consistent with the results of Zhu et al. 47 , which may be related to convection in the tropics.
ENSO exhibits seasonality 54 .Scatter plots of detrended water vapor anomalies and normalized IC2 (ENSO pattern) in winter, spring, summer, and autumn were generated to investigate the impact of ENSO on the upper tropospheric water vapor during different seasons (Fig. 9).The correlation coefficients between the upper tropospheric water vapor anomalies and normalized IC2 in winter, spring, summer, and autumn are 0.62, 0.74, 0.53, and 0.47, respectively.The results indicate that the correlation coefficient between upper tropospheric water  www.nature.com/scientificreports/vapor and ENSO pattern is higher during winter and spring, while it is lower during summer and autumn.This may be attributed to the seasonality of ENSO, as it peaks during DJF, but it takes months for the effects of ENSO to propagate to the upper tropospheric water vapor.The above results indicate that the influence extent of ENSO on upper tropospheric water vapor varies with the changing seasons.

Conclusion
Understanding the characteristic variations of upper tropospheric water vapor and identifying the controlling factors influencing them are crucial for comprehending global climate change.Although many studies have indicated the impact of ENSO and global warming on water vapor, the relationship of ENSO and global warming on the variations in upper tropospheric water vapor is not yet fully understood.ERA5 specific humidity data from the European Centre for Medium-Range Weather Forecasts (ECMWF) is a product of combining data assimilation with weather forecast modelling system, the accuracy of which relies on data assimilation methods and the stability and precision of observational systems.On this basis, we extracted the temporal and spatial patterns of ENSO and global warming from January 1959 to December 2021 by employing the strategy of ICA combined with the non-parametric optimal dimension selection to address the problem of the non-Gaussian distribution characteristics of climate time series and the uncertainty in the output dimension selection of ICA.
According to the results of global warming and ENSO patterns, the characteristics of water vapor variations in the upper troposphere and response to global warming and ENSO are revealed, providing crucial scientific insights for understanding global climate change.The main findings are as follows.
We analyzed the variations characteristics of water vapor in the upper troposphere and identified the factors affecting water vapor.There is an obvious seasonal variability of water vapor in the upper troposphere, which is usually the highest in summer and the lowest in winter.In addition, there is a clear upward trend in the time series of specific humidity anomalies, which may be associated with global warming according to the Clausius-Clapeyron equation.The water vapor in the upper troposphere has obvious latitude characteristics.The upper tropospheric water vapor is mainly concentrated in the tropics (30° N-30° S), and under the influence of the positive feedback mechanism of global warming, the longer the time scale, the greater the intensity and extent of the water vapor response in the upper troposphere.Most importantly, global warming and ENSO are crucial factors contributing to variations in upper tropospheric water vapor, with water vapor contributions due to global warming and ENSO patterns accounting for 86% of the global interannual variations in water vapor.In addition, the remaining 14% of the water vapor change may be associated with water vapor transport caused by wind, seasonal changes, and other factors.In particular, peaks in water vapor anomalies are associated with years of strong El Niño, while minimum values are associated with years of strong La Niña.The response of upper tropospheric water vapor to both global warming and ENSO takes a certain amount of time (i.e., it takes 11 months of propagation for the global warming pattern to have an impact on upper tropospheric water vapor, while ENSO takes 8 months of propagation).Global warming is the response to a steady, globally uniform forcing, thus taking a longer time compared to ENSO.
This study quantitatively analyzes the response relationship between variations in upper tropospheric water vapor and global warming.It is shown that not only does the global warming pattern (IC1) account for 53% of the global interannual variations in water vapor, but also a significant correlation with the upper tropospheric specific humidity anomalies is achieved.On the other hand, the spatial pattern of the global warming (IC1) is also very similar to the spatial pattern of the 63-year mean of upper tropospheric specific humidity.The above results indicate that global warming pattern (IC1) captures the majority of water vapor variations and is one of the primary factors driving variations in upper tropospheric water vapor.The correlation coefficient between the global warming pattern and the upper tropospheric water vapor anomalies is 0.87 and the global warming pattern also explains as much as 53% of the global interannual variations in water vapor, indicating that the upper tropospheric water vapor change is highly dependent on global warming.This may be related to the following reason.According to the Clausius-Clapeyron equation, atmospheric water vapor increases with rising temperatures 49 .Therefore, the atmospheric water vapor content will increase under the influence of global warming.Water vapor is also a greenhouse gas, and with the increase in greenhouse gases, the amount of longwave radiation reaching the Earth's surface also intensifies, further exacerbating global warming 23 .Consequently, this leads to a heightened correlation between global warming patterns and water vapor in the upper troposphere.Under the influence of global warming, there also will be significant changes in the water cycle at the global and regional scales, as well as an intensification of water vapor transport between the oceans and the land.Furthermore, both upper tropospheric water vapor anomalies and GMST achieve local peaks or minima during strong El Niño or strong La Niña years.This is related to the fact that strong ENSO events may have a significant impact on water vapor variations.Under the influence of the positive feedback mechanism of water vapor, it contributes to the variations of greenhouse gases, thereby influencing global warming.
The study also quantitatively analyzed the relationship between upper tropospheric water vapor variations and ENSO.It is shown that the ENSO pattern (IC2) accounts for 33% of the global interannual variations in water vapor and also has a strong correlation with the upper tropospheric specific humidity anomalies.Without considering the phase, the spatial pattern of the ENSO (IC2) is very close to the distribution of ENSO-year mean of upper tropospheric specific humidity during the DJF months.The above results suggest that ENSO captures some of the water vapor variations and is also one of the main factors contributing to water vapor variations in the upper troposphere.The results also show that the influence extent of ENSO on upper tropospheric water vapor varies with the changing seasons.During El Niño events, the distribution of global upper tropospheric water vapor is characterized by mainly positive anomalies, while during La Niña events, it is characterized by mainly negative anomalies.Additionally, relatively dry patches (blue patches) were observed near the warm pool during El Niño years, which is likely associated with tropical convection.

Figure 1 .
Figure 1.Processing flow chart of component extraction using the strategy of the optimal lowpass filtering combined with ICA.

Figure 2 .
Figure 2. Processing flow chart of data preprocessing (a-c) and the component extraction of the affecting factors (d, e).(a) Global mean specific humidity; (b) Deducting the annual cyclic and monthly variations; (c) Processing of time series with optimal lowpass filter based on (b); (d) GMST and IC1 extracted by the strategy of ICA combined with the non-parametric optimal dimension selection; (e) ONI and IC2 extracted by the strategy of ICA combined with the non-parametric optimal dimension selection.the years of strong El Niño and La Niña events (the intensity of ENSO events covered in this paper is provided by the website: https:// www.ggwea ther.com/ enso/ oni.htm) are highlighted in all the time series figure with light red and light blue backgrounds respectively.
depicts the latitude-time evolution of specific humidity anomalies data after filtering, from which we can also observe similar spatial characteristics of water vapor anomalies variations.In addition, we observe a notably strong positive response around the years 1972-73 (Strong El Niño), 1982-83 (Very Strong El Niño), 1987-88 (Strong El Niño), 1997-98 (Strong El Niño), 2010-11 (Strong La Niña), 2015-16 (Very Strong El Niño), and 2020-21 (Moderate La Niña).These years are almost all associated with strong ENSO events, indicating a significant impact of strong ENSO events on the upper tropospheric water vapor anomalies.

Figure 3 .
Figure 3. Latitude-time evolution of the global zonal mean specific humidity (Units: kg/kg).(a) Global specific humidity; (b) deducting the annual cyclic and monthly variation of the global specific humidity; (c) processing of time series with optimal lowpass filter based on (b).

Figure 4 .
Figure 4.The latitude distribution of the deseasonalized water vapor (WV) anomalies (black) and of the water vapor anomalies related to the global warming pattern (red), ENSO pattern (blue) and global warming plus ENSO pattern (green) derived from ICA.

Figure 5 .
Figure 5. Correlation analysis of the time series of IC1 and the upper tropospheric water vapor.(a) the time series of the upper tropospheric water vapor anomalies (red line) and IC1 (blue line) extracted by the strategy of ICA combined with the non-parametric optimal dimension selection; (b) the time series of GMST (red line) and IC1 (blue line).The time series in this paper are normalized to overcome the effect of amplitude uncertainty in separating components with ICA.

Figure 6 .
Figure 6.Analysis of the spatial relationship between global warming and the upper tropospheric water vapor.(a) The 63-year mean of upper tropospheric specific humidity at the 250 hPa pressure level (Units: kg/kg); (b) the spatial pattern of IC1 (Units: kg/kg).(c) The 63-year mean of the vertically integrated water vapor fluxes at global grid points (Units: mg/s kg).

Figure 7 .
Figure 7. Correlation analysis of the time series of IC2 and the upper tropospheric water vapor.(a) The time series of the upper tropospheric water vapor anomalies after removal of linear trends (red line) at the 250 hPa pressure level and IC2 (blue line) extracted by the strategy of ICA combined with the non-parametric optimal dimension selection; (b) the time series of ONI (red line) and IC2 (blue line).The time series in this paper are normalized to overcome the effect of amplitude uncertainty in separating components with ICA.

Figure 8 .
Figure 8. Analysis of the spatial relationship between ENSO and the upper tropospheric water vapor (Units: kg/kg).(a) The El Niño-year mean of upper tropospheric specific humidity at the DJF months at the 250 hPa pressure level; (b) the La Niña-year mean of upper tropospheric specific humidity at the DJF months at the 250 hPa pressure level; (c) the spatial pattern of IC2.

Figure 9 .
Figure 9. Scatter plots of water vapor anomalies versus normalized IC2 (ENSO pattern) in winter (DJF = December, January, and February), spring (MAM = March, April and May), summer (JJA = June, July and August), and autumn (SON = September, October, November) (Units: kg/kg).The values in the upper left corner are the correlation coefficients, where the p-values for correlation coefficients R in winter, spring, summer, and autumn are all less than 0.01.